Simple Microscopic Theory of Amontons' Laws for Static Friction 
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A microscopic theory for the ubiquitous phenomenon of static friction is presented. Interactions 
between two surfaces are modeled by an energy penalty that increases exponentially with the degree 
of surface overlap. The resulting static friction is proportional to load, in accordance with Amontons' 
laws. However the friction coefficient between bare surfaces vanishes as the area of individual 
contacts grows, except in the rare case of commensurate surfaces. An area independent friction 
coefficient is obtained for any surface geometry when an adsorbed layer of mobile atoms is introduced 
between the surfaces. The predictions from our simple analytic model are confirmed by atomistically 
detailed molecular dynamics simulations. 



C/3 



X3 
O 

> 

o 
o 
o 

-I— > 



I 

o 
o 



X 
S3 



Static friction Fs is the lateral force that must be ap- 
plied to initiate sliding of one object over another. Its 
presence implies that the objects have locked together 
into a local energy minimum that must be overcome by 
the external force. This threshold force is observed be- 
tween all the objects around us, yet its molecular origins 
have remained baffling. How do surfaces manage to lock 
together, and why does the friction obey Amontons' laws 
that Fs increases linearly with the load L pushing the ob- 
jects together, and is independent of the surface area [Q? 

Early theories of friction were based on the purely ge- 
ometric argument that friction is caused by interlocking 
of surface asperities The idea (Fig. |l|) is that the 

top surface must be lifted up a typical slope tana de- 
termined by roughness on the bottom surface. If there 
is no microscopic friction between the surfaces, then the 
minimum force to initiate sliding is — Liana. This 
result satisfies Amontons' laws with a constant coefficient 
of friction /ig = F^/L = tana. In 1737, Belidor obtained 
a typical experimental value of /is ~ 0.35 by modeling 
rough surfaces as spherical asperities arranged to form 
commensurate crystalline walls More recently, Is- 
raelachvili et al. have discussed a similar "cobblestone" 
model where the spherical asperities are atoms [|| . How- 
ever, asperities on real surfaces do not match as well as 
envisioned in these models or sketched in Fig. |l|. On av- 
erage, for every asperity or atom going up a ramp, there 
is another going down. One concludes that the mean fric- 
tion between rigid surfaces vanishes unless they happen 
to have the same periodicity and alignment. Detailed 
calculations show that elastic deformations are generally 
too small to alter this conclusion jj-j^]. 

In this Letter, we study a simple model of the static 
friction in individual contacts as a function of their area 
A. At fixed A, the static friction between bare surfaces 
is proportional to load. However, as implied by previous 
work [^-^, the prefactor vanishes for incommensurate 
surfaces. We find that F^/L also vanishes for disordered 
surfaces, scaling as as the area increases. These 

results for bare surfaces are inconsistent with Amontons' 
laws. Of course most surfaces around us are covered by 



a few angstroms of hydrocarbons, water and other small 
airborn molecules jsj. We show that introducing such 
mobile molecules into the interface between incommen- 
surate or disordered surfaces yields a value of /is = Fs/L 
that is independent of contact size and load. This natu- 
rally leads to Amontons' laws for any distribution of con- 
tacts. Our analytic predictions are then tested against 
previous |^ and new computer simulations. 

In our model, two surfaces "pay" an energy penalty 
V^vw that increases exponentially as the local separation 
between them decreases and they begin to overlap. This 
form of interaction is frequently used to model the Pauli 
repulsion between atoms, and contains the hard-sphere 
models of Belidor and Israelachvili et al. as a limiting 
case. Within the contact area A the surfaces are parallel 
to the xy plane. The bottom wall is held fixed, has zero 
mean height, and local height 5z\^{x). The mean height 
Zt and lateral position xt of the top wall are varied, and 
its local height is Zt + 5zt{x — Xt). Then T4fw can be 
written as an integral over A 



Kvw(a;t, ^t) 



(fa 



^-[zt+Szt(x-Xt)-Sz-b{x)]/i, 



(1) 



where ^ and e characterize the length and energy per unit 
area of the interaction. The dependence of T^w on zt can 
be pulled out of the integral in Eq. |l]. The integral then 
depends only on xt and can be expressed as an effective 
shift in the mean wall separation by lS.z{xt)'. 



K;w(^t, Zt) = sA e 



-{zt+Az{xt)}/(, 
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The lateral and normal forces on a static contact are 
equal to the derivatives of V„w with respect to Xt and 
Zt, respectively. Due to the simple form of Eq. ^ almost 
all factors in the ratio of friction to normal load cancel, 
and the ratio is independent of load |^ . The static fric- 
tion coefficient corresponds to the maximum of this ratio, 
taken along the direction of the applied force: 



fis = max 



(3) 
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Physically, dAz{xt)/dxt represents an effective ramp 
that the top surface must climb to escape the local poten- 
tial energy minimum. It is analogous to tana in Fig. 
but arises from many molecular-scale interactions. Note 
that the factorization that leads to Eq. |^ relies on the 
exponential form of the local energy penalty. However 
our numerical tests with other interaction potentials and 
the simulations with Lennard- Jones potentials described 
below show that the result is more general. 

The effect of the area of contact on the friction coeffi- 
cient fj,s is most easily seen by making a cumulant expan- 
sion of Eq. (|l|). To leading order, Az is proportional to 
the mean-squared difference between Szt and Sz^. Only 
the contribution from the cross product contributes to 
the variation of Az with respect to xt , yielding 



max — ikx5z\,{k)5zl{k)e 

k 



ik-xt 
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where 6zh,t sue the 2D Fourier transforms of (5zb,t- Ex- 
plicit evaluation of Eq. ^ for various trial surfaces, and 
the simulations described below, show that higher order 
terms in Sz do not alter any of the following predictions. 

Eq. (^) implies several relations between the geome- 
try of the two surfaces and the friction between them: 
(i) If the two surfaces are crystalline and commensurate, 
then Szt and Szh share a common periodicity. The cor- 
responding Bragg peaks in their Fourier transforms lead 
to an area independent fig- (ii) The magnitude of fig is 
largest for two identical surfaces, where all Bragg peaks 
contribute. For commensurate surfaces that are not iden- 
tical, /is decreases exponentially with the length of the 
common period. For example, if there are p lattice con- 
stants of the bottom surface for every q lattice constants 
of the top surface, then the contribution to /ig comes 
from the qth Bragg peak of the bottom surface and the 
Pth Bragg peak of the top surface. The Fourier content 
of atomically rough surfaces [|lo| drops at least exponen- 
tially with |fc|, producing a corresponding decrease in /Xg- 
(iii) As a consequence of (ii), fj,s vanishes completely for 
infinite incommensurate contacts. Contributions from 
the circumference of finite contacts yield a rapidly van- 
ishing contribution to /ig (iv) For two disordered 
but smooth surfaces Szt{k) and Szi,{k) have rings of dif- 
fuse scattering that overlap. There will be many k that 
contribute to fig, but the phase of each contribution will 
be random. For interfaces with contact-area independent 
corrugation {Sz^-^), one can immediately conclude from 
the law of large numbers that /ig c>c A~'^/'^. Self-similar 
and curved interfaces need a more careful treatment pT| . 

Two geometry independent predictions can also be 
drawn from Eqs. (j^) and (^: (v) /ig does not depend on 
the interaction strength e. (vi) Allowing for elastic defor- 
mation of the surfaces typically reduces /Xg, because the 
roughness decreases as the surfaces become more compli- 
ant 0. 



Predictions (i) to (iii) agree with previous analytic and 
simulation studies of the force between commensurate 
and incommensurate surfaces p[-|q|. (v) is in agreement 
with recent computer simulations Iq] , where doubling the 
strength of Lennard- Jones interactions produced almost 
no change in /ig. (iv) and (vi) are new results that are 
tested by simulations below. First we consider the impli- 
cations of mobile atoms between the surfaces within our 
analytic model. 

Fig. 1^ illustrates how a mobile layer of atoms between 
two surfaces can lead to a finite /ig independent of the 
geometry of the two surfaces. The mobile atoms are able 
to move to positions where they simultaneously match 
the geometry of both top and bottom surfaces. This has 
the effect of augmenting the height of the bottom surface 
in a way that matches the undulations of the top surface, 
and gives an area independent contribution to Eq. ^. 

In order to incorporate mobile atoms into our theory 
quantitatively, we presume that wherever a mobile atom 
sits in the interface, the effective distance between top 
and bottom wall is reduced by the diameter d of the 
mobile atom. This effectively incorporates the mobile 
atoms into one of the walls. The position of a mobile 
atom in the interface is defined by a vector Xi so that the 
dimensionless areal density of the mobile atoms can be 
written as p{x) = {ird^ /A) 5{x — Xi). For the sake of 
simplicity, we assume that the mobile layer screens the 
direct wall-wall interaction and that interactions between 
the mobile atoms can be ignored. The latter assumption 
is consistent with observations that /ig is insensitive 
to the density of mobile atoms. The validity of both as- 
sumptions was verified through simulations for the model 
system introduced below. 

With this model, the total energy of the system is given 
by the indirect wall-wall interaction mediated through 
the mobile atoms. This takes the form 

Kvw = d'x p{x) e-{-t+5.t(S-?t)-5.t(5)-d}/c 
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The force on the top wall Ft can again be calculated by 
taking the gradient of Kvw with respect to the center of 
mass position of the top wall: (a?t,zt)- It is possible to 
decompose Ft into forces Ft.i that stem from individual 
atoms. For each such force, a linear relationship between 
the X component of the force Ft^ix and the "local" load 
Ftiz can be established that is similar to Eq. (0): 



Ft_ 



_d_ 

dxt 



{Szt{xi ~ Xt) ~ 6zh{xi)) 



F, 



(6) 



As long as the partial derivative has a well-defined av- 
erage value, the total friction will rise linearly with the 
total load. 
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At this point, we turn back to a more qualitative dis- 
cussion. As long as the temperature is small compared 
to the energy barrier for diffusion between the surfaces, 
atoms will sit near local energy minima. The simulations 
described below show that these local minima usually 
correspond to the (++) configurations identified in Fig. 
^, where the energy with respect to both surfaces is con- 
cave. Lateral displacements of the walls will ultimately 
limit the volume accessible to atoms and increase 

V-ujw The atoms have created interlocking asperities that 
resist sliding much like those in Fig. |^. One can show 
more rigorously from Eq. || that if all atoms sit at (++) 
positions, the response to a lateral displacement is an op- 
posing force that is linear in the normal load L and inde- 
pendent of the area of contact. Indeed, this "Amontons' 
law for elastic pinning" with an elastic restoring force 
linear in L has been observed experimentally [p^ . Thus, 
we may conclude that mobile atoms in the interface lead 
to Amontons' laws with a non-vanishing friction coeffi- 
cient. The pinning of two surfaces through a "between- 
sorbed" layer can also be interpreted in terms of a gen- 
eralized Frenkel-Kontorova or Tomlinson model Q. In 
these models, incommensurate layers exhibit static fric- 
tion when one of the surfaces is so compliant that it can 
conform to the other surface. The mobile atoms act like 
an elastic sheet with nearly infinite compliance and thus 
allow locking to occur. 

We have tested the above predictions for disordered 
walls using molecular dynamics simulations of the model 
described in detail in Refs. ||^ and |Q. The two walls 
contain discrete atoms tied to their equilibrium sites by 
springs of stiffness n. Periodic boundary conditions are 
applied in the plane of the walls. Mobile molecules be- 
tween the surfaces contain n spherical monomers bound 
into a chain. All monomers and wall atoms also interact 
with a Lennard- Jones potential. For the results presented 
below, n = 1 or 6, and the energy and length scales (e 
and cr, respectively) of the Lennard-Jones potential are 
the same for monomers and wall atoms. Other parame- 
ters produced equivalent results. Disordered walls were 
made by quenching bulk fluid states and then removing 
all atoms above or below some height. 

Fig. H shows representative results for the scaling of 
friction with contact area. The top wall was coupled 
to a constant normal load and a slowly increasing lateral 
force. The static friction was determined from the thresh- 
old force needed to initiate lateral motion of the top wall. 
For all cases considered, the coefficient of friction be- 
tween bare disordered surfaces vanished as A~^^^ (solid 
line) in agreement with prediction (iv). In contrast, when 
molecules were inserted into the interface, the static fric- 
tion rapidly approached a constant as A was increased. 
It is interesting to note that the bare and disordered 
surfaces have equal friction at an area corresponding to 
about 1000 atoms. This is slightly larger than contacts in 
atomic force microscopy, but much smaller than the mi- 



cron size contacts between typical surfaces. In agreement 
with prediction (vi), decreasing k to the smallest reason- 
able value [|j (k = 400 ea^^) reduced the friction between 
bare walls considerably (crosses). Much smaller reduc- 
tions (~20%) were observed when mobile atoms were 
present between elastic surfaces. In both cases the re- 
ductions reflect the ability for atoms to move vertically 
to minimize the steric overlap with the opposing surface. 
When atomic displacements were constrained to the hor- 
izontal plane, we observed no reduction in Fg. 

The simulations also allowed us to test the assumption 
that monomers in (H — \-) configurations (Fig. |^) dominate 
the pinning. Fig. ^ shows results for a quarter monolayer 
of spherical molecules {n — 1) as the force was ramped up 
to a value slightly above Fg over a time iicvorsc and then 
returned four times as rapidly to zero. About 90% of 
mobile atoms sit at (-1- -I-) positions until the wall begins 
to slide (i/ircvcrso ~ 0.95), and they contribute an even 
larger fraction of Fg. Atoms that sit at convex positions 
relative to both walls (--), were only seen during sliding. 
Even in this dynamic state, almost 70% of the atoms are 
in (+ +) sites and they continue to provide most of the 
lateral force. A detailed analysis of these runs shows that 
the mechanisms of kinetic and static friction are closely 
related in this model. At any instant in time, most of 
the mobile atoms are in local energy minima and resist 
lateral motion. When these sites become unstable (+-) 
or (- -) sites, the atoms pop rapidly to a new site. Energy 
is dissipated during these rapid pops and flows into the 
walls as heat. 

In conclusion, we have presented a simple analytic 
model for the molecular origins of friction. Although the 
model assumes exponential repulsion between surfaces, 
we presented simulations with more realistic potentials 
that show that the predictions are more generally ap- 
plicable. Specifically the model provides a microscopic 
foundation for Amontons' law i^iatcrai = fJ-L, but shows 
that /i vanishes for bare incommensurate or disordered 
surfaces as the size of contacts increases. Introducing mo- 
bile atoms between the surfaces yields Amontons' laws: 
fi is independent of surface area and load for any contact 
geometry. The simulations also reveal deep connections 
between static and kinetic friction. Of course, experi- 
mental systems contain additional features that have not 
been treated. These include long-range elastic and plas- 
tic deformations of the walls, generation of wear debris, 
and chemical reactions [|. All of these may be important 
in particular systems. 
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FIG. 1. Sketch of two surfaces with interlocking asperities. 
The top surface experiences a normal load L and a lateral 
force F, which attempts to pull the top surface up the slope 
tana. The bottom wall is fixed. 




FIG. 2. Sketch of two rigid, incommensurate surfaces that 
are separated by a submonolayer of mobile atoms. Circles in- 
dicate positions where the gap between upper and lower wall 
is a local maximum. The larger the gap, the smaller the en- 
ergy penalty for occupation by a mobile atom. Full circles 
indicate (-1- -I-) positions where atoms fit into concave regions 
of both surfaces. The open circle indicates a less favorable 
(H — ) site where one wall is concave and the other is convex. 
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FIG. 3. Static friction coefficient as a function of con- 
tact area A. The friction coefficient for rigid or elastic 
(k = 400e(T~^) bare walls drops as A~^^^ (solid line). In- 
serting enough mobile molecules (n = 6) to make either a 
quarter or two monolayers yields nearly the same area inde- 
pendent value of Us- Statistical error bars are comparable to 
the symbol size. 
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FIG. 4. (a) Fraction of atoms in (+ +), (+-), and (--) sites 
(Fig. 2) and (b) their contribution to the restoring friction 
force as a function of time. For these runs treverse was about 
300 in Lennard- Jones units. 
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